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ABSTRACT 

We report the discovery of a 1 niHz optical quasi-periodic oscillation (QPO) in the 
candidate ultracompact low-mass X-ray binary 4U 06144-091. The ultra-low frequency 
QPO has no X-ray counterpart in contemporaneous RXTE/PGA data and is likely a 
signature of structure in the accretion disk. The QPO can be reasonably fitted with a 
single sine wave but with a phase jump part way through the observation, indicating 
that it is not coherent. We also identify a 48 m in mo dulation, approximately consistent 
with the suggested orbital period of lO'Brienl (|2005[ ) and Shahbaz et al., (,2008.) . If this 
is indeed orbital, it supports an identification of 4U 0614-1-091 as an ultra-compact 
source. 
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1 INTRODUCTION 

Low-mass X-ray binaries (LMXBs) can be divided into 
black hole candidates and neutron star systems. Those 
with orbital periods of hours to days are expected to have 
hydrogen-rich mass donors. A few ultra-compact X-ray bi- 
naries (UCXBs) containing hydrogen-deficient and possi- 
bly degenerate donors can evolve to extremely small bi- 
nary sep arations, with orbital periods a s short as a few 
minutes (|Nelson. Rappaport. &: Jos3ll986l '). Although these 
1 systems were initially assumed to be rare, many can- 
didates have now been fou nd ( Ncl emans fc JonkeJ l2006l : 

lin't Zand. Jonker. fc Markw ardt 2007). 

The X-ray source 4U 0614 +091 was first optically iden- 
tified bv iMurdin et all (|l974l ') with a faint (~ 18 mag), 
blue, var iable star (V1055 Ori) located in the Galac- 
tic plane (|van Paradiis fc van der "Kii3ll994l ). Several argu- 
ments recently concl uded 4U 0614-1-091 is a UCXB host - 
ing a neutron star. lJuett. Psaltis. fc Chakrabartvl (|200lh 
found an enhanced neon to oxygen ratio in 4U 0614-1-091 
similar to that seen in ot her UCXBs. T h rough sen- 
sitive optical spectroscopy, iNelemans et all (|2004l ) and 

INelemans. Jonker. fc Steeghd (|2006l ') found significant car- 
bon and oxygen emission lines, but no hydrogen or helium. 
They argued 4U 061 4-1-091 contai n s a ca rbon-oxygen accre- 
tion disk. Recently, iMadei et all (|2010l ^ reported the dis- 
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covery of a broad emission feature at ~ 0.7 keV, which 
is attribute of O VIII Lya emission, in the spectra of 
4U 0614-1-09. This f eature has been seen so far in two systems 
iMadei et al.ll201ll ) that are receiving oxygen-rich material 
from the CO or ONe white dwarf. It is an endorsement of 
an oxygen-rich mass donor star. Thermonuclear typ e I X-ray 
bursts were observed by OSO-8 JSwank et all 19781'). a nd as - 
suming they were Eddington limited. iBrandt et al.l (|l992l ) 
deduced a distance to 4U 0614- 1-091 probably < 3kpc. More 
recently iKuulkers et al.l (|2010h have re-examined archival 
bursts from several satellites and derived a distance of 
3.2 kpc. These authors argue that the bursts require more 
helium in the accreted matter than would be consistent with 
the optical spectra leaving the nature of the mass donor still 
in question. 

A straightforward resolution of this puzzle has been 
elusi ve due to the difficu lty in identifying an orbital pe- 
riod. iMachin et al] j 19901 ) found a 10-d photometric mod- 
ulation which might be caused by precession of an accre- 
tion disc. They then scaled the precession period of Her X- 
1 to 4U 0614-1-091 to infer the orbital period. If the ra- 
tio of the disc radius to the neutron star Roche lobe ra- 
dius, p is assumed to be ~ 1, Porb ~ 6.3 — 9.4 hr. If 
p ~ 0.5 then Porb ~ 2.4 — 3.3 hr. Either way, this argu- 
ment suggests a non-UCXB system with a hydrogen-rich 
donor. However, the faintness of its X-ray emission in spite 
of a nearby distance suggest 4U 0614-1-091 has a very low 
acc retion rate more consiste nt with an orbital period < 1 
hr ijPelove fc BildstenI [200^ ) . Recently [O'BrienI ||2005|) re- 
ported a tentative 50 min optical period and Shahbaz et al.l 
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l|2008l ) also claimed a strong modulation of 51.3 mins. A 
spectroscopi c modulation on a period close to th is may also 
be present ijNelemans. Jonker. fc Steeehd 120061 ) . These re- 
cent results support the case that 4U 0614+091 is a UCXB 
system. 

The orbital period is not the only periodicity expected. 
Superhumps at close to the orbital period can occur in 
short period (and usually small mass ratio) LMXBs due 
to the beat between the orbital period and a slowly pre- 
cessing accretion disk (|Haswell et akllioOll ) . Millihertz quasi- 
periodic oscillations (QPOs) also occur, but are rarely ob- 
served at optical wavelengths. The exception is 4U 1626- 
67 which is a LMXB source containing a 7.66 s X-ray pul- 
sar in an ultra-compact 42min orbit. 130 mHz ('^ 7.7s) 
X-ray pulsations are seen, a 48 mHz (~ 20.8s) QPO is 
produced by an interaction between the pulsar's magneto- 
sphere and the accreti on disk, and a 1 niHz (~ 16.7 min) 
UV QPO were found l|Chakrabartv et al.ll200lh . Only the 
1 mHz QPO was not detected in simultaneous X-ray obser- 
vations. IChakrabartv et aP (j200lh suggested that the 1 mHz 
oscillations are due to warping of the inner accretion disk. 

We present here new optical photometry of 
4U 0614-1-091 revealing a new mHz QPO and further 
evidence for a ~ 50 min period. In Section [2] we describe 
our data. The lightcurves are analysed in Section |3] and 
compared with contemporaneous X-ray observations in 
Section |4l Finally in Section [5] we discuss the implications 
of a ~ 50 min orbital period and the possible origin of the 
QPO. 



Table 1. Optical observing log for 4U 0614+091. 



Name 


Observation Date 


Duration 
(sec) 


Exp Time 

(sec) 


Night 1 


12/08/2007 


27330 


10 


Night 2 


12/09/2007 


27450 


10 


Night 6 


12/13/2007 


26540 


10 



.5 1.0 
1 0.9 



Night 1 



5.0x10' 1.0x10* 1.5x10* 2.0x10* 2.5x10* 



Night 2 



5.0x10' 1.0x10* 1.5x10* 2.0x10* 2.5x10* 

Night 6 



5.0x10' 1.0x10* 1.5x10' 2.0x10* 2.5x10' 

Time (sec) 



Figure 1. Light curves of 4U 0614+091 during the three nights 
(top panel: night 1; middle panel: night 2; bottom panel: night 
6). All fluxes are measured relative to a comparison star and 
normalized to unit mean flux. 



2 OBSERVATIONS AND DATA REDUCTION 

We performed fast optical photometry of 4U 0614+091 on 
three nights in December 2007 with the 2.1m (82") Otto 
Struve Telescope at McDonald Observatory, Texas (Ta- 
ble [1]). Approximately 7 hours of data were obtained on 
each night. For convenience we will refer to the nights as 
Night 1, Night 2, or Night 6 for the re mainder of this work. 
We u sed the Argos CCD Photometer (|Mukadam fc Natheil 
l2005h . This is a CCD camera designed for fast photometry 
that was originally designed for pulsating white dwarfs, but 
has proven ideal for interacting binaries. It is possible to ob- 
tain exposures as short as 1 second with negligible dead-time 
between exposures. Hence it is ideally suited to relatively 
short period phenomena, such as UCXB orbital periods and 
mHz QPOs. 

As listed in Table [1] a total of 8132 CCD image frames 
were taken with an exposure time of 10 sec per frame. Tim- 
ing was synced to GPS 1 s ticks and checked against multiple 
NTP servers. All images used a broad BVR filter to maxi- 
mize the source count rate. 

Custom IDL software written for Argos data was used 
to subtract dark current and bias, and apply flat fleld cor- 
rections to the data before performing aperture photometry 
of 4U 0614+091 and a brighter comparison star. All results 
reported here are based on differential photometry relative 
to the comparison star, which showed no signs of variability. 



3 OPTICAL ANALYSIS AND RESULTS 

3.1 Lightcurves 

Figure [1] shows light curves of 4U 0614+091 for our three 
observations. While much of the variability appears to be 
aperiodic flickering, an oscillation with a ~ 1000 sec period 
is prominent throughout night 2 and possibly in shorter seg- 
ments on the other two nights. 

3.2 Periodicity Analysis 

To better quantify the periodicities present, fast Fourier 
transforms (FFTs) were computed to derive power spec- 
tra. These power spectra are shown in Figure [2] Two ma- 
jor peaks are recurrent across multiple nights. A lower fre- 
quency at 0.34 mHz mainly appeared on Nights 1 and 6, 
and is coincident with earlier clai ms of a tentative Porb 
(|O'Brienll2005l : IShahbaz et al.ll2008l ). The higher frequency 
signal at 1.03 mHz is strong only on Night 2, but a weaker 
signal at a similar frequency is also detectable on the last 
night. The two peaks were fitted with Gaussians in order 
to measure characteristic frequencies. Fitting parameters 
are recorded in Tabled The 1.03 mHz peak corresponds 
to a 16.2 ± 0.1 min period and the 0.34 mHz peak to a 
49.0 ±0.1 min period. 

3.3 Multi-Lorentzian Analysis 

Where a mixture of periodic, quasi-periodic, and aperiodic 
variability is present, as may be the case in these data, a 
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Table 3. Best Fit Parameters for multiple Lorentzian Model. Average of Feature 1 considers Night 1 and Night 6. Average of Feature 2 
includes Night 2 and Night 6. 



Number 


White Noise 


VLF Noise 




Feature 1 






r 

(10-3) 


^max 

mHz 


r 

(10-3) 


mHz 


^max 

mHz 


r 

(10-3) 


Night 1 
Night 2 
Night 6 


2.66±0.06 
4.84±0.08 
5.56±0.13 


0.25±0.13 
0.02±0.04 
0.03±0.02 


0.93±0.02 
0.12±0.23 
O.OliO.Ol 


0.35±0.01 
0.32±0.08 
0.34±0.02 


0.37±0.01 
0.46±0.03 
0.37±0.01 


0.62±0.03 
0.41±0.03 
0.42±0.02 


Average frequency 








0.35 ±0.01 


a 




Average period 








48 ±2 min 






Number 




Feature 2 






Feature 3 






mHz 


mHz 


r 

(10-3) 


mHz 


J^max 

mHz 


(10-3) 


Night 1 
Night 2 
Night 6 


1.17±0.01 
l.OSiO.Ol 
0.99±0.03 


1.18±0.01 
1.03±0.01 
l.OlitO.Ol 


0.10±0.03 
0.43±0.06 
O.lliO.Ol 


3.37±0.15 
2.94±0.10 
2.64±0.50 


3.64±0.15 
3.46±0.04 
3.76±0.20 


0.09±0.02 
0.10 ±0.01 
0.14±0.01 


Average frequency 


1.03 ±0.01^ 






3.46 ±0.04 


c 




Average period 


16.2 ±0.2 min 






4.8 ±0.1 min 





Average of night 1 and night 6 only. 
Average of night 2 and night 6 only. 
For broad features, Umax is considered. 




Table 2. Results of fitting Gaussians to linearly scaled power 
spectra of 4U 0614+091. 
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Figure 2. Linearly scaled power spectra of 4U 0614+091 for three 
observations (top panel; night 1; middle panel: night 2; bottom 
panel: night 6). 



superposition of multiple Lorentzians has become widely 
used for power spe c tral a n alysis in LMXBs (Nowak 200 ^; 
van _Straaten et"all I2OO2I : iBelloni. Psaltis. fc van der KlisI 
While these techniques were developed to describe 
X-ray data, there is no reason that they should not also be 
applicable to optical data too. The primary advantage of a 



Name 




Peak Freq 
(mHz) 


FHWM 
(mHz) 


Period 
(min) 


Night 1 


Peak 1 


0.34 + 0.01 




48 ±2 


Night 2 


Peak 2 


1.03 ±0.01 


0.10 ±0.02 


16.1 ±0.2 


Night 6 


Peak 3 
Peak 4 


0.34 ±0.12 
1.07 ±0.12 


0.09 ± 0.03 
0.19 ±0.17 


48 ±3 
15 ± 2 


Low Freq Avg 
High Freq Avg 


Peak 1&3 
Peak 2&4 


0.34 ±0.01 
1.03 ±0.01 




49.0 ±0.1 
16.2 ±0.1 



Lorentzian description of the variability components is that 
all quasi-periodic oscillations (QPOs) and noise components 
can be fitted in the same way; it is not necessary t o decide 
the na ture of each component in advance. Following lNowal3 
(|2000l ). an individual Lorentzian component is given by. 



1 



(1) 



TV A2 + (l/- 1^0)2 

where r is the integrated fractional r.m.s. of each Lorentzian, 
and A is its half width at the half maximum (HWHM). 

We show in Figure |3] the power s pectra and th e fitting 
functions in the (fPu) representation l|Nowakl|200 0l). where 
the power spectral density is multiplied by its Fourier fre- 
quency. It is helpful to associate a characteristic frequency 
with each Lorentzian component. We follow X-ray conven- 
tion and use Vmax as the characteristic frequency for broad 
features. This is the maximum in the i/P^ representation 
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and defined by, 

I^max = = 1^0 ^1+^ (2) 

see iBelloni. Psaltis. &: van der Klij l|2002h . where Q is the 
quahty factor, defined as Q = t'o/2A. For narrow features. 

Previous research using this technique on different 
sources was done in a higher frequency range, 0.01- 
100 Ifz, typical of frequ en cies seen in X-ray data, e.g. 
Ivan Straaten et al. 

I (I2OO2'), Bclloni, Psalt is. &: van der KlisI 
( 20021 ). These authors have found correlations between char- 
acteristic frequencies and break frequencies of low and high 
frequency noise. Our optical results show similar complexity, 
with multiple components required, but in a lower frequency 
range than typically studied in X-ray observations, exploit- 
ing the longer uninterrupted time-series obtainable from the 
ground. 

We found that to obtain a satisfactory fit, four 
Lorentzian components plus statistical white noise were 
needed. We began with the two frequencies already isolated, 
but found that we could not obtain a fit without excess 
power remaining at both higher and lower frequencies. One 
zero-centered Lorentzian component was added for very low 
frequency noise (VLF noise); this is a band-limited noise 
component. The last Lorentzian component is fitting the 
higher frequency peak revealed only in the i/P^, represen- 
tation as residual power around 3mHz. The characteristic 
frequencies and integrated noise are listed in Table [31 The 
quoted errors in z^o, i^max and r are la confidence limits. 

As clearly seen in Figure O the three peak features can 
all be identified on each night in this representation. The 
first two of them are relatively narrow features. Because no 
significant differences are present between uo and Vmax, we 
take Umax as the frequency of the feature. Feature 1 with an 
average of 0.35 mHz is most significant on Night 1 and Night 
6. Hence, we ignore results from Night 2 for this period. The 
second sharp feature with average of 1.03 mHz is especially 
strong on Night 2 but is also present on Night 6. Results from 
Night 1 are ignored for Feature 2. Generally, the average 
results for these features are very consistent with results 
from Gaussian fitting of the linearly scaled power spectra. 
Feature 3 is quite broad. The average of Umax from all three 
nights is calculated. It is centered at 3.46 mHz with about a 
4% uncertainty. 



3.4 Phase stability of the 1 mHz QPO 

Unlike the other features, the 1 mHz QPO can be clearly 
seen in the lightcurves, allowing direct study in the temporal 
domain. In order to better understand the 1 mHz QPO, we 
examined the stability of its phase. To do this we performed 
a sliding sine wave fit to the Night 2 data. We used a 1000 sec 
window shifting with small step sizes of 100 sec and fitted 
a sine wave fixed at the 16.2 min period, but with its phase 
and amplitude as a free parameter. The change in phase 
with time is shown in Figure|4l An unexpected phase change 
occurs around 9500 sec. The phase is basically separated into 
two segments, differing by 0.34 in phase. As a check, we also 
examined individual segments of the light curve by eye and 
confirmed that the phase change can be seen at this time. 




0.01 0.10 1.00 10.00 100.00 

Frequency (mHz) 
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Figure 3. Power spectra and fit functions in the power spectral 
density times frequency representation for three observations (top 
panel: night 1; middle panel: night 2; bottom panel: night 6). The 
different dotted lines mark the individual Lorentzian components 
of the fit. The dark line displays the final model. 



This confirms that the signal is quasi-periodic in nature and 
not a true period. 
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Figure 4. Phase Plot of Night 2 Period 1.03 mHz using 1000 sec 
as window size and step is 100 sec. The error bars are measured 
for each 1000 sec window. 

4 COMPARISON WITH X-RAY DATA 

Our observations lie outside of the fre quency range com- 
monly studied in X-rays. For example, Ivan Straaten et"aD 
l|2002h terminate the power spectra densities (PSDs) of 
4U 0614-1-091 around 10 mHz. Some of their PSDs do re- 
veal a very low frequency noise component overlapping with 
our frequency range. 

We have searched for an X-ray mHz QPO with contem- 
poraneous X-ray observations from the Rossi X-Ray Timing 
Explorer (RXTE). We chose the Proportional Counter Ar- 
ray (PCA) B band data, which is concentrated in a low 
energy band, 2-5 KeV and examined data from dates near- 
est to our observation, from Dec 13-15 2007. Lightcurves 
were produced in the usual way using HEASoft version 6.8. 
Lightcurves were extracted in 16 sec bins, with a maximum 
duration of 51.2 min. 

We used the lightcurves to generate low-frequency 
power-spectra, and show these in Fig. [5] There is no ev- 
idence for an X-ray oscillation coincident with our 1 mHz 
optical feature, although this may be a consequence of the 
non-simultane ity of the data. This is sim ilar to the situation 
in 4U 1626-67 ijChakrabartv et al.ll200ll ). This is a 42 min or- 
bital period system showing a strong 15 min UV and optical 
QPO without an X-ray counterpart. 



5 DISCUSSION 

Including the new results here, there are now four 
independent observations supporting an optic al pe- 
riodicit y of app roximately 50 min ('O'Brien '2005': 
[Nelcm ans. Jonker. fc Stecghs 2006; Shahbaz ct al. 2008 and 
our results). Besides confirming the nature of 4U 0614-1-091 
as an ultracompact source, identification of the 50 min 
period as orbital would have important implications for the 
nature of the companion star. We should keep in mind, 
however, that a periodicity that is detected repeatedly is not 
necessarily an orbital period. The possible presence of it in 
spectroscopic data would sug gest an orbital origin, but this 
detec tion was very marginal ( Nelemans. Jonker. fc Steeghd 
l2006l ). 



1 10 
Frequency (mHz) 



Figure 5. Power-spectrum of RXTE PCA B band data from 
Dec 13-15 2007. 



iNelemans et al-l l|2004 ) suggested that the companion 
was a CO white dwarf based on optical detections of carbon 
and oxygen lines, but an absence of hydrogen and helium. 
INelemans et al.l (|2010l ') reconsidered this argument based on 
chemical abundance predictions from new models and favour 
either a hybrid white dwarf or a very evolved helium star 
donor. For a 50 min orbital period, these two possibilities 
become observationally similar. 

An additional consideration is that the presence of type 
I X-ray bursts s uggests the presence of helium in the ac- 
creted material (|Kuulkers et al.l l2010t ). Both helium star 
and hybrid white dwarf models could include mass-transfer 
phases where CO-enriched material with a non-negligible 
helium fraction is transferred, although the upper limits 
on the helium fraction implied by the non-detection of he- 
lium lines seem in contra diction to the helium fraction re- 
quired for X-ray bursts (|Kuulkers et al.ll2010t ). The pres- 
ence of apparent helium bursts in hel ium deficient UCXBs 
appe ars to be a recurrent problem (|juett fc Chakrabartvl 
120031 ). These authors proposed a generic solution to this 
problem in which CO rich material is accreted but then 
spallation reactions near the surface of the neutron star 
break these elements up in to hydrogen and helium which 
can then be burnt back to heavie r elem ents during bursts 
(jBildsten. Salpeter. fc Wassermanlll992l '). We note that this 
mechanism will not work in pure CO material as spalla- 
tion requires elements of quite different atomic number to 
be present. If some helium is present in the accreted mate- 
rial, however, then this mechanism could potentially work 
to increase the helium fraction above that seen in the disk 
spectra. 

A third possibility for a 50 min binary would be that 
the donor star is an evolved main-sequence star that has 
lost most of its hydrogen envelope. In this case, however, 
we would expect CNO processed abundances with nitrogen 
more abundant than c arbon and oxygen, in contradiction to 
optical spectroscopy INe lemans et al.ll2004l ). so this possi- 
bility can probably be discounted. 

A major difficulty with identifying the donor as a white 
dwarf or evolved helium star comes in the expected X- 
ray luminosity of the system. X-ray observations indicate 
an average luminosity of 3.6 x lO^^ergs"^. In a neutron 
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star system with no event horizon the X-ray luminosity 
should be a reasonable indicator of mass transfer rate, and 
in any case a lower-limit on it. One uncertainty remaining 
is that 4U 0614+091 does eject jets (Migliari et_al. 2006). 
A more comprehensive study, however, suggests that only 
10 p ercent of the accretio n luminosity is being carried by the 
jets (|Migliari et al.|[2010l ). This constraint, together with the 
observed X-ray luminosity, implies a mass transfer rate of 
3 X 10"^" Mq yr~^. This is much larger (by an order of mag- 
nitude) than expected for the hybrid white dwarf or evolved 
heli um star models at 50 min and evolving to longer peri- 
ods (jNelemans et al.ll2010l ). The inferred mass transfer rate 
would be about right for an evolved main-sequence star, but 
as discussed above this appears inconsistent with the chem- 
ical abundances. 

It is certainly possible that the 50 min periodicity is not 
the orbital period, as none of the measurements demonstrate 
that the deduced modulations have to be orbital in origin. 
A white dwarf or helium star scenario would give the right 
mass transfer rate f or a shorter period be low 30mins (as 
earlier suggested bv^Nelemans e t al]|2004l '). In this case, our 
mHz QPO would be a candidate for th e orbital period, an d 
might then be attributed a super hump (|Haswell et al.ll200ir ). 
It is probably not the orbital period which should be strictly 
coherent, as we have seen phase shifts. Even in the case of a 
superhump, however, it is hard to see how a phase shift such 
as this could arise since it would involve a readjustment of 
the accretion disk on close to its dynamical timescale. If the 
orbital period is indeed below 30 min, we then also have to 
explain the origin of the recurrent 50 min timescale. 

A more natural explanation for a large amplitude mHz 
oscillation seen in the optical but not X-rays, and that is 
able to undergo large phase shifts in a single cycle is that it 
originates further into the accretion disk but in regions too 
cool to produce X-rays. Optical light in 4U 0614+091 ap- 
pears to be dominated by the accretion disk (|Migliari et al.l 
I2OO6I . I2OIOI ) so participation by the disk is probably neces- 
sary to obtain a 10 percent amplitude. The frequency and 
strong amplitude of our mHz o scillation strikingly resem- 
bles that seen in 4U 1626-67 IChakrabartv et al.l [2OO1I ). 
That system too shows a CO spectrum and is identified 
as like ly containing a hybri d white dwarf or helium star 
donor ((Nelemans et al.ll2010l) . The orbital period of 42 mm 
is a close match to that proposed for 4U 0614+091, but 
much more secure since it is seen in pulse-timing residu- 
als. For this orbital period the mass transfer rate inferred 
for a white dwarf donor is also about an order of magni- 
tude greater than th at expected from gravitational radiation 
IChakrabartvl (Il998 '). a problem similar to that encountered 
if we attribute a 50 min orbital period to 4U 0614+091. 

IChakrabartv et al] l|200ll ) suggested that the oscilla- 
tion in 4U 1626-67 arises from a purely geometric ef- 
fect, the tilt angle of the normal to the local disk surface 
varies with azimuth. The reprocessed emission would be ex- 
pected t o be sensitive to warpin g material along the line 
of sight. IShirakawa fc Tail (|2002l l suggested the magnetic 
torque from the pulsar interacting with the accretion disk 
may produce the warping effect. The similarity of behav- 
ior in 4U 0614+091, however, appears at odds with this as 
there is no evidence that this system also harbors a neutron 
star with a sig nificant m agnetic field. A possible solution 
is proposed bv iLail (|2003l ). In the absence of a magnetized 



compact object, it is suggested that magnetic fields thread- 
ing the disk itself could cause magnetically driven warping. 
This idea was proposed in the context of systems with a 
magnetically-driven outflow. Multiwavelength observations 
suggest the presence of a jet-like outflow is 4U 0614+091 
(|Migl iari et al. I I2OO6I . I2OIOD and so this may be a natural 
system in which to invoke this effect. 

As a final remark, we note that the mHz oscillation 
is not seen alone, but as one of three frequencies that are 
recurring from night to night. While one of these may be 
orbital, it is possible that none are. We are thus not looking 
for a mechanism that can give rise to a single oscillation, 
but to a more complex multi-periodic behavior, rather rem- 
iniscent of that see at m uch higher frequencies in X-rays 
(|van Straaten et al.ll200j ). This would suggest a more com- 
plex multi-modal warping of the accretion disk as a possible 
mechanism. We have focused on the 1 mHz QPO here be- 
cause it can be seen clearly in the lightcurves, and so studied 
in the time-domain. The 3.5 mHz QPO is certainly also of 
interest but with less information, and no obvious timescale 
to associate it with, it is harder to say anything about it. 
More observations will be needed to properly disentangle 
the behavior seen, for example by looking for correlations 
between changes in the oscillation frequencies to test if they 
are linked or independent. 



6 CONCLUSIONS 

We have identified three persistent, distinct, periodic 
or quasi-periodic oscillations in the optical flux from 
4U 0614+091. The 0.35 mHz (48 min) p eriod was suggested 



to be an orbital period by'C^'Brion (200 5|) and Shahbaz et al 
( 20081) , and supported by Nelemans. Jonker. fc Steeghs 
(120061 ). If it is indeed orbital, it would confirm 4U 0614+091 



as an ultra compact binary, (jjuett. Psaltis. fc Chakrabartvl 
I2OOII : iNelem ans et al. 200i). 

The origin of the 1 mHz QPOs now seen in both 
4U 0614+091 and 4U 1626-67 remains unknown, ff they 
are even related to each other. It likely originates in the ac- 
cretion disk and could arise as the precession period of inner 
disk regions that are not hot enough to produce X-ray emis- 
sion. The less prominent 3 mHz QPO may be related to this, 
or may be distinct. 
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